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A brief summary of several topics in the study of gravitational many body problem is given. The 
discussion covers both static backgrounds (applicable to astrophysical systems) as well as clustering 
in an expanding background (relevant for cosmology). 
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INTRODUCTION 

The statistical mechanics of systems dominated by 
gravity is of interest both from the theoretical and "prac- 
tical" perspectives. Theoretically, this field has close con- 
nections with areas of condensed matter physics, fluid 
mechanics, renormalization group, etc. From the practi- 
cal point of view, the ideas find application in different ar- 
eas of astrophysics and cosmology, especially in the study 
of globular clusters, galaxies and gravitational clustering 
in the expanding universe. [For a review of statistical me- 
chanics of gravitating systems in static background, see 
0; textbook descriptions arc in 2|; gravitational cluster- 
ing in cosmology is reviewed in and in the textbooks 
for a sample of different approaches see (5(1 and the 
references cited therein.] 



GRAVITATIONAL CLUSTERING IN STATIC 
BACKGROUND 

To construct the statistical description of a system of 
A'' self gravitating point particles, one should begin with 
the construction of the micro canonical ensemble describ- 
ing such a system. If g{E) is the volume of the constant 
energy surface H{pi^qi) — E, then the entropy and the 
temperature of the system will be S{E) = lng{E) and 
T{E) = (3{E)-^ = (dS/dE)-^. (The finiteness of g re- 
quires the system to be confined to a finite volume in 
space for any system). 

Systems for which a description based on canonical 
ensemble is possible, the Laplace transform of g{E) with 
respect to a variable /3 will give the partition function 
Z[P). Gravitating systems of interest in astrophysics, 
however, cannot be described by a canonical ensemble 
Qi IS' Virial theorem holds for such systems and we 
have {2K + U) =0, where K and U are the total ki- 
netic and potential energies of the system. This leads to 
E = K + U = —K; since the temperature of the system is 
proportional to the total kinetic energy, the specific heat 
will be negative: Cy = {dE/dT)v oc [dE/dK) < 0. On 
the other hand, the specific heat of any system described 
by a canonical ensemble Cy = /3^((Ai?)^) will be posi- 
tive definite. Thus, one cannot describe self gravitating 
systems of the kind we are interested in by canonical en- 



semble. 

One can, however, attempt to find the equilibrium con- 
figuration for self gravitating systems by maximising the 
entropy S{E) or the phase volume g{E). For a system of 
point particles, there is again no global maximum for en- 
tropy 0] , • If we move a small number of these particles 
arbitrarily close to each other, the potential energy of in- 
teraction of a pair of these particles, —Gmim2/ri2, will 
become arbitrarily high as ri2 — > 0. Transferring some 
of this energy to the rest of the particles, we can increase 
their kinetic energy without limit. This will clearly in- 
crease the phase volume occupied by the system (in the 
momentum space) without bound. This argument can be 
made more formal by dividing the original system into a 
small, compact core, and a large diffuse halo and allowing 
the core to collapse and transfer the energy to the halo. 

The absence of the global maximum for entropy - as 
argued above - depends on the lack of small scale cutoff. 
If we assume, instead, that each particle has a radius a, 
there will be an upper bound on the amount of energy 
that can be made available to the rest of the system. Fur- 
ther, no real system is completely isolated and to obtain 
a truly isolated system, we need to confine the system 
inside a spherical region of radius R with, say, reflecting 
wall. 

The two cut-offs a and R will make the upper bound 
on the entropy finite, but even with the two cut-offs the 
formation of a compact core and a diffuse halo will still 
occur, since this is the direction of increasing entropy. 
Particles in the hot diffuse component will permeate the 
entire spherical cavity, bouncing off the walls and hav- 
ing a kinetic energy which is significantly larger than the 
potential energy. The compact core will exist as a grav- 
itationally bound system with very little kinetic energy. 
A formal way of understanding this phenomena is based 
on the virial theorem [ij : 



2T + U ^ 3PV + $ 



(1) 



for a system with a short distance cut-off confined to a 
sphere of volume V, where P is the pressure on the walls 
and is the correction to the potential energy arising 
from the short distance cut-off. This equation can be 
satisfied in essentially three different ways. If T and U 
are significantly higher than 3PV and then we have 
2T + U ^ which describes a self gravitating systems 
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in standard virial equilibrium, but not in the state of 
maximum entropy. If T ^ U and SPV ^ one can 
have 2T Ki 3PV which describes an ideal gas with no po- 
tential energy confined to a container of volume V; this 
will describe the hot diffuse component at late times. If 
T <^ U and iPV <C then one can have C/ de- 
scribing the compact potential energy dominated core at 
late times. Such an asymptotic state with two distinct 
phases is quite different from what would have been ex- 
pected for systems with only short range interaction. If 
the gravitating system is put in a heat bath and the tem- 
perature is varied, a sudden phase transition occurs at a 
critical temperature, leading to the formation of the two 
phases 0, 0- 

There are, however, configurations which are local ex- 
trema of entropy, which are not global maxima. Intu- 
itively, one would have expected the distribution of mat- 
ter in such configuration to be described by a Boltzmann 
distribution, with the /9(x) oc exp[— /30(x)], where (j) is the 
gravitational potential related to the density p by Poisson 
equation. This configuration, called isothermal sphere, 
has a density profile p oc x~'^ asymptotically. Isother- 
mal spheres with total energy E and mass M , however, 
cannot exist Q if {RE/GM'^) < -0.335. Even when 
(RE/GAP) > -0.335, the isothermal solution need not 
be stable. The stability of this solution can be investi- 
gated by studying the second variation of the entropy. 
Such a detailed analysis shows that the following re- 
sults are true: (i) Systems with {RE/G]VP) < -0.335 
cannot evolve into isothermal spheres. Entropy has no 
extremum for such systems 

Q, 0- (ii) Systems with 
{{RE /GAP) > -0.335) and (p(0) > 709 p{R)) can exist 
in a meta-stable (saddle point state) isothermal sphere 
configuration. Here p{0) and p{R) denote the densities 
at the center and edge respectively. The entropy extrema 
exist but they are not local maxima, (iii) Systems with 
{{RE/GM^) > -0.335) and (p(0) < 709 p{R)) can form 
isothermal spheres, which are local maximum of entropy. 
These are striking peculiarities in the case of SMGS and 
seem to find application in the physics of globular clus- 
ters. 



GRAVITATIONAL CLUSTERING IN AN 
EXPANDING BACKGROUND 

There is considerable amount of observational evidence 
to suggest that one of the dominant energy densities in 
the universe is contributed by self gravitating (nearly) 
point particles. The smooth average energy density of 
these particles drive the expansion of the universe while 
any small deviation from the homogeneous energy den- 
sity will cluster gravitationally. It is often enough (and 
necessary) to use a statistical description and relate dif- 
ferent statistical indicators (like the power spectra, nth 
order correlation functions, ....) of the resulting den- 



sity distribution to the statistical parameters (usually the 
power spectrum) of the initial distribution. 

The relevant scales at which gravitational clustering is 
nonlinear are less than about 10 Mpc, while the expan- 
sion of the universe has a characteristic scale of about 
4000 Mpc 4]. Hence, nonlinear gravitational clustering 
in an expanding universe can be adequately described by 
Newtonian gravity by introducing a proper coordinate for 
the i— th particle r^, related to the comoving coordinate 
Xi, by Ti — a{t)x.i where a{t) is the expansion factor. 
The Newtonian dynamics works with the proper coordi- 
nates which can be translated to the behaviour of the 
comoving coordinate by this rescaling. 

If x(t, q) is the position of a particle at time t with its 
initial position being q, then equations for gravitational 
clustering in an expanding universe, in the Newtonian 
limit, can be summarised as [3- 

.. 2d . 1 , 

x-K— x = — ^Vx0; (2) 

where x = x(t, q), g(k, p) = (fc/2)2 + p'^ - 2(k.p/fc)2 
and 0k (i) is the Fourier transform of the gravitational 
potential (/)(t, x) due to perturbed density. 

Equation Q is exact but involves x(t, q) on the right 
hand side and hence, cannot be considered as closed. To- 
gether, the two equations form a closed set but solving 
them exactly is an impossible task. It is, however, possi- 
ble to use this equation with several well motivated ap- 
proximations to obtain information about the system. 
I shall briefly mention a few of them. 

Consider flrst the effect of a bunch of particles, in a 
virialized cluster, on the rest of the system. This is de- 
scribed, to the lowest order, by just the monopole mo- 
ment of the cluster - which can be taken into account 
by replacing the cluster by a single particle at the cen- 
tre of mass having appropriate mass. Such a replacement 
should not affect the evolution at scales much bigger than 
the cluster size. At first sight, one may wonder how this 
feature ( "renormalizability of gravity") is taken care of 
in equation 0. Inside a galaxy cluster, for example, 
the velocities x can be quite high and one might wonder 
whether this could infiuence the evolution of 0k at all 
scales. This does not happen and, to the lowest order, 
the contribution from virialized bound clusters cancel , 
[ifij l in the two terms in the right hand side of ©• 

In this limit, 0k is constant in time and the Pois- 
son equation — fc^0k = 47rGp;,a^5k oc {Sy^/a) in the 
matter dominated universe with a{t) oc t^^^,pi, oc 
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implies that the density contrast has the growing so- 
lution (5k(i) — [a{t) / a{ti)]6it,{ti) . The power spectrum 
P(k, t) — (|5k(^)|^) and the correlation function ^(x, i) 
[which is the Fourier transform of the power spectrum] 
both grow as a^{t). This allows us to fix the evolution of 
clustering at sufficiently large scales uniquely. The clus- 
tering at these scales, which is well described by linear 
theory, grows as a^. 

There is, however, an important caveat to this claim. 
While ignoring the right hand side of © one is compar- 
ing its contribution at any wave number k to the contri- 
bution in linear theory. If at the relevant wavenumber, 
the contribution from linear evolution is negligibly small, 
then the only contribution will come from the terms on 
the right hand side and, of course, we cannot ignore it in 
this case. This contribution will scale as k^B?, where R 
is the typical scale of virialized systems and will lead to 
a development of (5k oa k^,P{k) oc fc^ at small k. Thus, 
if the large scales have too little power intrinsically (i.e., 
if n > 4 ), then the long wavelength power will soon 
be dominated by the "fc^ - tail" of the short wavelength 
power arising from the nonlinear clustering. This is an 
interesting and curious result which is characteristic of 
gravitational clustering. 



NONLINEAR SCALING RELATIONS 

As to be expected, cosmological expansion completely 
changes the asymptotic nature of the problem. The prob- 
lem has now become time dependent and it will be point- 
less to look for "equilibrium solutions" . 

There are three key theoretical questions which are 
of considerable interest in this area which I will briefly 
summarise: 

• If the initial power spectrum is sharply peaked in 
a narrow band of wavelengths, how does the evo- 
lution transfer the power to other scales? (This is, 
in some sense, analogous to determining the Green 
function for the gravitational clustering except that 
superposition will not work in the nonlinear con- 
text.) 

• Do the virialized structures formed in an expand- 
ing universe due to gravitational clustering have 
any invariant properties? Can their structure be 
understood from first principles? 

• Does the gravitational clustering at late stages wipe 
out the memory of initial conditions or does the late 
stage evolution depend on the initial power spec- 
trum of fluctuations? 

To make any progress with these questions we need a 
robust prescription which will relate statistical indicators 
like the two-point correlation function in the nonlinear 



regime to the initial power spectrum. Fortunately, this 
problem has been solved jl^ to a large extent and hence, 
one can use this as a basis for attacking these questions. 

The nonlinear mean correlation function can be ex- 
pressed in terms of the linear mean correlation function 
by the relation: 



D 



11.7^(0,0 



Dh/2 



(for f < 1) 

(for 1< <e < 125) (4) 
(for 125 < 



where / — x[l + ^{a,x)]^^^ , = 2,3 is the dimension 
of space and h is a, constant. [The results of numerical 
simulation in 2D, suggests that /i = 3/4 asymptotically. 
We will discuss the 3D results in more detail below] . The 
numerical values are for Z? = 3. 

One could use this to examine whether the power spec- 
trum (or the correlation function) has a universal shape 
at late times, independent of initial power spectrum. 
This is indeed true [13] if the initial spectrum was sharply 
peaked. In this case, at length scales smaller than the ini- 
tial scale at which the power is injected, the two point 
correlation function has a universal asymptotic shape of 
f (a, x) oc a^x^^{L + x)^^, where L is the length scale at 
which ^ « 200. This can be understood as follows: 

In the quasi-linear phase, regions of high density con- 
trast will undergo collapse and in the nonlinear phase 
more and more virialized systems will get formed. We re- 
call that, in the study of finite gravitating systems made 
of point particles and interacting via Newtonian grav- 
ity, isothermal spheres play an important role and are 
the local maxima of entropy. Hence, dynamical evolu- 
tion drives the system towards an (l/cc^) profile. Since, 
one expects similar considerations to hold at small scales, 
during the late stages of evolution of the universe, we may 
hope that isothermal spheres with (1/x^) profile may still 
play a role in the late stages of evolution of clustering 
in an expanding background. However, while converting 
the density profile to correlation function, we need to dis- 
tinguish between two cases. In the quasi-linear regime, 
dominated by the collapse of high density peaks, the den- 
sity profile around any peak will scale as the correlation 
function and we will have ^ oc (l/x^). On the other hand, 
in the nonlinear end, we will be probing the structure in- 
side a single halo and ^(x) will vary as {p{x + y) p{y)) . 
If p oc [xl"*^, then ^ oc with 7 = 2e — 3. This gives 

^ oc (l/x) for e = 2. Thus, if isothermal spheres are 
the generic contributors, then we expect the correlation 
function to vary as (1/x) and nonlinear scales, steepening 
to (1/x^) at intermediate scales. Further, since isother- 
mal spheres are local maxima of entropy, a configuration 
like this should remain undistorted for a long duration. 
This argument suggests that a f which goes as (1/x) at 
small scales and (l/x^) at intermediate scales is likely to 
grow approximately as at all scales. At scales bigger 
than the scale at which power was originally injected, the 
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spectrum develops a fc** tail for reasons described before. 
This is confirmed by simulations for sharply peaked ini- 
tial spectra. But if the initial spectrum is not sharply 
peaked, each band of power evolves by this rule and the 
final result is a lot messier. 

The second question one could ask, concerns the den- 
sity profiles of individual virialized halos. If the density 
field p(a,x) at late stages can be expressed as a super- 
position of several halos, each with some density profile 
/(x) then the i-th halo centred at x^ will contribute a 
density /{k — Ki,a) at the location x. The power spec- 
trum for the density contrast, S{a,x) = (p/pb — l), will be 
P{k) = |/(fc)pPc(fc), where Pc(k,a) denotes the power 
spectrum of the distribution of centers of the halos. If 
the correlation function ^ oc x the correlation func- 
tion of the centres ^ oc x^'^" and the individual profiles 
are of the form f{x) oc a;"*^, then this relation translates 
toe = 3+(l/2)(7-7,). 

At very nonlinear scales, the centres of the virialized 
clusters will coincide with the deep minima of the grav- 
itational potential. Hence, the power spectrum of the 
centres will be proportional to the power spectrum of 
the gravitational potential P,p{k) oc fc""^ if P{k) oc fc". 
Since the correlation functions vary as when the 

power spectra vary as fc", it follows that 7 = 7c — 4. 
Substituting into the above relation, we find that e = 1 
at the extreme nonlinear scales. On the other hand, in 
the quasi-linear regime, reasonably large density regions 
will act as cluster centres and hence, one would expect 
Pc(fc) and P{k) to scale in a similar fashion. This leads 
to 7 ~ 7c, giving e w 3. So we would expect the halo 
profile to vary as x~^ at small scales steepening to x~^ 
at large scales. A simple interpolation for such a density 
profile will be 



Such a profile, usually called NFW profile ^^Ij is often 
used in cosmology. The argument given above, however, 
is very tentative and it is difficult to obtain |(SJ) from a 
more rigorous theoretical analysis. 

In fact, it is possible to reach different conclusions re- 
garding the asymptotic evolution of the system from dif- 
ferent physical assumptions Q| . The NSR in Q for 3-D 
with constant h leads to the asymptotic correlation func- 
tion 

F/ \ ^ --y 3/i(n + 3) 

2 + h[n + 3) 

for an initial spectrum which is scale-free power law with 
index n. If we assume that the evolution gets frozen 
in proper coordinates at highly nonlinear scales then it 
is easy to show that h = 1. If this assumption (called 
stable clustering) is valid, then the late time behaviour 
of f (a, x) is strongly dependent on the initial conditions 



and lO shows that £,{a,x) at nonlinear scales will be as, 

^(a,x) oc a^x"^^; (f» 200). (7) 

In other words the two (apparently reasonable) require- 
ments: (i) validity of stable clustering at highly non- 
linear scales and (ii) the independence of late time be- 
haviour from initial conditions, are mutually exclusive. 
[At present, there exists some evidence from simulations 
|15l | that this process, called stable clustering, does not 
occur in the a oc t^/"^ cosmological model; but this result 
is not definitive]. 

In the very nonlinear limit, the correlation function 
probes the interiors of individual halos and we have e = 
(l/2)(3-|-7). [This corresponds to Pc = constant, 7c = 3 
in the previous discussion.] If 7 depends on n so will e 
and the individual halos will remember the initial power 
spectrum. 

We can obtain a 7 which is independent of initial power 
law index provided h satisfies the condition h{n + 3) — c, 
a constant. In this case, the halo profile will be given by 
e = 3(c-|-l)/(c-|-2). Note that we are now demanding the 
asymptotic value of h to explicitly depend on the initial 
conditions though the spatial dependence of f (a, x) does 
not. As an example of the power of such a — seemingly 
simple — analysis, note the following: Since c > 0, it 
follows that e > (3/2); invariant profiles with shallower 
indices (for e.g with e = 1) discussed above are not con- 
sistent with the evolution described above. One requires 
very high resolution simulations to verify the condition 
h{n + 3) = c and the current results are inconclusive. 
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